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We present a concise review of the recent theoretical progress concerning the standard model calcu- 
lation of the inclusive radiative B — > Xsj decay. Particular attention is thereby devoted to the cal- 
culations of the next-to-next-to-leading order fixed-order 0{a'l) contributions, non-local 0{asA/mi,) 
power corrections, and logarithmic-enhanced 0{a.g) cut-effects to the decay rate. The current status 
of new physics calculations of the inclusive fe 57 mode is also briefly summarized. 
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I. 



INTRODUCTION 



As a flavor-changing-neutral-current process the in- 
clusive radiative i?-nieson {B = B° or B^) de- 
cay is Cabibbo-Kobayashi-Maskawa- (CKM) and loop- 
suppressed within the standard model (SM) and thus 
very sensitive to new physics (NP) effects. In order to 
exploit the full potential of i? ^ ^s7 iu constraining 
the parameter space of beyond the SM physics both the 
measurements and the SM prediction should be known 
as precisely as possible. 

The present experimental world average (WA) which 
includes the latest measurements by CLEO, Belle, and 
BaBar [1] is performed by the Heavy Flavor Averaging 
Group [2] and reads for a photon energy cut of Ej > Ecut 
with Ecut — 1.6 GeV in the B-meson rest-frame 

B{B ^ X,7)exp = (3.55 ± 0.24l° °|; ± 0.03) x lO""* . (1) 

The total error of the WA is below 8% and consists 
of i) a combined statistical and systematic error, ii) a 
systematic uncertainty due to the extrapolation from 
^-cut = [1.8, 2.0] GeV to the reference value, and in) a 
systematic error due to the subtraction of the B Xd^ 
event fraction. At the end of the B-factory era the final 
accuracy of the averaged experimental value is expected 
to be around 5%. 



II. BASIC PROPERTIES OF B 



Xs-r 



The b sj transition is dominated by perturba- 
tive QCD effects which replace the power-like Glashow- 
Iliopoulos-Maiani (GIM) suppression present in the elec- 
troweak (EW) vertex by a logarithmic one. This mild 
suppression of the QCD corrected amplitude reduces the 
sensitivity of the process to high scale physics, but en- 
hances the B ^ Xsj branching ratio (BR) with respect 



t Based on invited talks given at the third meeting of the workshop 
"Flavour in the era of the LHC", CERN, Geneva, March 15-17, 
2006 and at XLIInd Rencontres de Moriond, QCD and Hadronic 
Interactions, La Thuile, Italy, March 17-24, 2007. 



to the purely EW prediction by a factor of around three. 
The logarithmic GIM cancellation originates from the 
non-conservation of the tensor current which is generated 
at the EW scale by loop diagrams involving ly-boson and 
top quarks exchange. The associated large logarithms 
L ~ hiMiy/mh have to be resummed at each order in 
Qfs, using techniques of the renormalization group (RG) 
improved perturbation theory. Factoring out the Fermi 
constant Gp, the b —f s"f amplitude receives corrections 
of 0«L") at leading order (LO), of 0«L"-i) at next- 
to-leading order (NLO), and of 0{a^L"-~'^) at next-to- 
next-to-leading order (NNLO) in QGD. 

A suitable framework to achieve the necessary resum- 
mation is the construction of an effective theory with 
five active quarks, photons and gluons by integrating out 
the top quark and the EW bosons. Including terms of 
dimension up to six in the local operator product expan- 
sion (OPE) the relevant effective Lagrangian at a scale fi 
reads 



QCDxQED 



4Gf 

7f 



fc=i 



(2) 



Here the first term is the conventional QCD and QED 
Lagrangian for the light SM particles. In the second term 
Vij denotes the elements of the CKM matrix and Ck (/^) 
are the Wilson coefficients of the corresponding operators 
Qk built out of the hght fields. 

The operators and the numerical values of their Wilson 
coefficients at ^ rrib are given by 



?i.2-(sr.c)(cr^6), 

h-6 ^ mb)j:,mq) , 

h = m^SLO^^bnE^, , 
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Ci,2{rnb) ~ 1 , 
1^3-6(^6)1 < 0.07, 
C7(mfc) ~ -0.3, 
C8(mb) ~ -0.15, 



(3) 



where F and F', entering both the current-current oper- 
ators Qi^2 and the QCD penguin operators Q3-6, stand 
for various products of Dirac and color matrices [3]. In 
the dipole operator Q7 (Qs); e [g) is the electromag- 
netic (strong) coupling constant, q^^R are the chiral quark 
fields, Ff^u is the electromagnetic (gluonic) field 

strength tensor, and T° are the color generators. 

After including LO QCD effects the dominant contri- 
bution to the partonic decay rate stems from charm quark 
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loops that amount to ^ 158% of the total b ^ sj decay 
amplitude. The top contribution is compared to the one 
from the charm quark with ^ —60% less than half as 
big and has the opposite sign. Diagrams involving up 
quarks are suppressed by small CKM factors and lead at 
the amplitude level to an effect of a mere ^2%. 

All perturbative calculations of 6 — > 57 involve three 
steps: i) evaluation of the initial conditions Ck{fJ-w) of 
the Wilson coefficients at the matching scale /i^- ~ M,^ 
by requiring equality of Green's functions in the full 
and the effective theory up to leading order in (external 
momenta) /M„,, ii) calculation of the anomalous dimen- 
sion matrix (ADM) that determines the mixing and RG 
evolution of Cfe(^) from down to the B-meson scale 
fib ~ nib, and iii) determination of the on-shell matrix 
elements of the various operators at fib ~ mi,. Due to 
the inclusive character of the B Xgj mode and the 
heaviness of the bottom quark, mf, ^ A Aqcd, non- 
perturbative effects arise in the last step only as small 
corrections to the partonic decay rate. 



III. 



THEORETICAL PROGRESS IN B 



Xs-r 



At the NNLO level, the dipole and the four-quark op- 
erator matching involves three and two loops, respec- 
tively. Renormalization constants up to four loops must 
be found for b sj and b sg diagrams with four-quark 
operator insertions, while three-loop mixing is sufficient 
in the remaining cases. Two-loop matrix elements of the 
dipole and three-loop matrix elements of the four-quark 
operators must be evaluated in the last step. 

The necessary two- and three-loop matching was per- 
formed in [4] and [5] . The mixing at three loops was de- 
termined in [6] and at four loops in [7] . The two- loop ma- 
trix element of the photonic dipole operator together with 
the corresponding bremsstrahlung was found in [8] and 
subsequently confirmed in [9]. These calculations have 
been very recently extended to include the full charm 
quark mass dependence [10]. The three- loop matrix el- 
ements of the current-current operators were derived in 
[11] within the so-called large- /3o approximation. A cal- 
culation that goes beyond this approximation employs an 
interpolation in the charm quark mass [12]. The effect 
of still unknown NNLO contributions is believed to be 
smaller than the uncertainty that has been estimated af- 
ter incorporating the above corrections into the SM cal- 
culation [12, 13]. To dispel possible doubts about the 
correctness of this assumption, calculations of the miss- 
ing pieces are being pursued. 

The most impressive bit of the various NNLO calcula- 
tions is the one of the four-loop ADM that describes the 
0{al) mixing of the four-quark into the dipole operators 
[7] . It has involved the computation of more than 20000 
four-loop diagrams and required a mere computing time 
of several months on around 100 CPU's. 

Another crucial part of the NNLO calculation is the 
interpolation in the charm quark mass performed in [12]. 
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FIG. 1: Renormalization scale dependences of B{B Xs'y)sM 
at LO (dotted lines), NLO (dashed lines), and NNLO (solid lines) 
in QCD. The plots show from top to bottom the dependence on 
the matching scale fiw, the iJ-meson scale fit, and the charm 
quark mass renormalization scale fic- 



The three- loop 0{a^) matrix elements of the current- 
current operators contain the charm quark, and the 
NNLO calculation of these matrix elements is essential 
to reduce the overall theoretical uncertainty of the SM 
calculation. In fact, the largest part of the theoretical 
uncertainty in the NLO analysis of the BR is related to 
the definition of the mass of the charm quark [14] that 
enters the 0{as) matrix elements {s'y\Qi^2\b). The lat- 
ter matrix elements are non-vanishing at two loops only 
and the scale at which nic should be normalized is there- 
fore undetermined at NLO. Since varying tUc between 
iTic{mc) ^ 1.25 GeV and mdmb) ^ 0.85 GeV leads to a 
shift in the NLO BR of more than 10% this issue is not 
an academic one. 

Finding the complete NNLO correction to {sj\Qi^2\b) 
is a formidable task, since it involves the evaluation of 
hundreds of three-loop on-shell vertex diagrams that are 
presently not even known in the case rric = 0. The ap- 
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proximation made in [12] is based on the observation that 
at the physical point rric ~ 0.25 the large rric 'S> mi, 
asymptotic form of the exact 0{as) [3] and large-/?o 
0(0^^0) [11] result matches the small nic ^ Wb expan- 
sion rather well. This feature prompted the analytic cal- 
culation of the leading term in the rric 2> mt expansion 
of the three-loop diagrams, and to use the obtained in- 
formation to perform a interpolation to smaller values of 
rric assuming the O{a'^(3o) part to be a good approxima- 
tion of the full 0{a'^) result for vanishing charm quark 
mass. The uncertainty related to this procedure has been 
assessed in [12] by employing three ansatze with different 
boundary conditions at nic = 0. A complete calculation 
of the 0{al) corrections to {sj\Qi,2\b) in the latter limit 
or, if possible, for nic ^ 0.25 nib, would resolve this am- 
biguity and should therefore be attempted. 

Combining the aforementioned results it was possible 
to obtain the first theoretical estimate of the total BR 
of B — > Xgj at NNLO. For the reference value E'cut = 
1.6 GeV the result of the improved SM evaluation is given 
by [12, 13] 

6(5^ Xs7)sM = (3.15 ±0.23) X 10"'', (4) 

where the uncertainties from hadronic power corrections 
(5%), parametric dependences (3%), higher-order pertur- 
bative effects (3%), and the interpolation in the charm 
quark mass (3%) have been added in quadrature to ob- 
tain the total error. 

The reduction of the rcnormalization scale depen- 
dences at NNLO is clearly seen in Fig. 1. The most pro- 
nounced effect occurs in the case of the charm quark mass 
rcnormalization scale fic that was the main source of un- 
certainty at NLO. The current uncertainty of 3% due to 
higher-order effects is estimated from the variation of the 
NNLO curves. The central value in Eq. (4) corresponds 
to the choice fJ,w,b,c = (160, 2.5, 1.5) GeV. More details 
on the phenomenological analysis including the list of in- 
put parameters can be found in [12]. 

It is well-known that the OPE for B Xg^ has cer- 
tain limitations which stem from the fact that the pho- 
ton has a partonic substructure. In particular, the local 
expansion does not apply to contributions from opera- 
tors other than Qj, in which the photon couples to light 
quarks [15, 16]. While the presence of non-local power 
corrections was thus foreseen such terms have been stud- 
ied until recently only in the case of the (QsiQs) hi- 
terference [15]. In [17] the analysis of non-perturbative 
effects that go beyond the local OPE have been extended 
to the enhanced non-local terms emerging from {Q7, Qs) 
insertions. The found correction scales like 0{asA/mb) 
and its effect on the BR was estimated using the vacuum 
insertion approximation to be —[0.3,3.0]%. A measure- 
ment of the flavor asymmetry between X^^ and 
B^ could help to sustain this numerical esti- 
mate [17]. Potentially as or maybe even more impor- 
tant than the latter correction are those arising from the 
{Qi,2tQ7) interference. Naive dimensional analysis sug- 
gests that some non-perturbative corrections to them also 



scale like ©(a^A/mb). Since at the moment there is not 
even an estimate of those corrections, a non-perturbative 
uncertainty of 5% has been assigned to the result in 
Eq. (4). This error is the dominant theoretical uncer- 
tainty at present and thought to include all known [17] 
and unknown 0{as A/nib) terms. Calculating the precise 
impact of the enhanced non-local power corrections may 
remain notoriously difficult given the limited control over 
non-perturbative effects on the light cone. 

A further complication in the calculation of B Xg^ 
arises from the fact that all measurements impose strin- 
gent cuts on the photon energy to suppress the back- 
ground from other i3-meson decay processes. Restricting 
to be close to the physical endpoint iSmax = ™b/2, 
leads to a breakdown of the local OPE. which can be 
cured by resummation of an infinite set of leading-twist 
terms into a non-perturbative shape function [18]. A de- 
tailed knowledge of the shape function and other sublead- 
ing effects is required to extrapolate the measurements to 
a region where the conventional OPE can be trusted. 

The transition from the shape function to the OPE 
region can be described by a multi-scale OPE (MSOPE) 
[19]. In addition to the hard scale jih ^ ™& ~ 5 GeV, this 
expansion involves a hard-coUinear scale ^hc ~ y/nib/S. ~ 
2.5 GeV corresponding to the typical hadronic invari- 
ant mass of the final state Xj, and a soft scale Hs ^ 
A ~ 1.5 GeV related to the width A/2 = mf,/2 - E^nt 
of the energy window in which the photon spectrum is 
measured. In the MSOPE framework, the perturbative 
tail of the spectrum receives calculable corrections at all 
three scales, and may be subject to large perturbative 
corrections due to the presence of terms proportional to 
as(VWA) 0.27 and a^(A) ~ 0.36. 

A systematic MSOPE analysis of the {Q7, Q7) interfer- 
ence at NNLO has been performed in [20]. Besides the 
hard matching corrections, it involves the two-loop log- 
arithmic and constant terms of the jet [19, 22] and soft 
function [23] . The three- loop ADM of the shape function 
remains unknown and is not included. The MSOPE re- 
sult can be combined with the fixed-order prediction by 
computing the fraction of events 1 — T that lies in the 
range Ecnt = [1-0, 1.6] GeV. The analysis [20] yields 

1-T = 0.07t°;°^p^^^ ± 0.02hadr ± 0.02pars , (5) 

where the individual errors are perturbative, hadronic, 
and parametric. The quoted value is almost twice as large 
as the NNLO estimate 1-T= 0.04±0.01pcrt obtained in 
fixed-order perturbation theory [12, 13, 21] and plagued 
by a significant additional theoretical error related to low- 
scale perturbative corrections. These large residual scale 
uncertainties indicate a slow convergence of the MSOPE 
series expansion in the tail region of the photon energy 
spectrum. Given that A is always larger than 1.4 GeV 
and thus fully in the perturbative regime this feature is 
unexpected. 

Additional theoretical information on the shape of the 
photon energy spectrum can be obtained from the uni- 
versality of soft and coUincar gluon radiation. Such an 



4 



Model 


Accuracy 


Effect 


Bound 


THDM type II 


NLO [27, 28] 


f 


Af± > 295 GeV (95% CL) [13] 


MFV MSSM 


NLO [28-33] 


t 


— 


LR 


NLO [28] 


t 


— 


general MSSM 


LO [34] 




\{5i3)LL\<4xlO-\ |(5|3)flfl| < 8 X 10-\ 
l(^l3)Lfl|<6xlO-^ \{5i,)RL\< 2x10-^ 


mUED 


LO [36] 




l/R > 600 GeV (95% CL) [41] 


RS 


LO [37] 




Mkk > 2.4 TeV 


LH 


LO [38] 


T 




LHT 


LO [39] 


I 





TABLE I: Theoretical accuracy, effect on B{B Xs^) relative to the SM prediction, and if applicable, constraint on the parameter 
space following from B Xs'y in popular NP scenarios. Arrows pointing upward (downward) indicate that the NP effects interfere 
constructively (destructively) with the SM b s'y amplitude. Single (double) arrows specify whether the maximal possible shift is 
smaller (larger) than the theoretical uncertainty of the SM expectation. See text for details. 



approach can be used to predict large logarithms of the 
form ln(iS'max — i^cut)- These computations have also 
achieved NNLO accuracy [24] and incorporate Sudakov 
and renormalon resummation via dressed gluon exponen- 
tiation (DGE) [24, 25]. The present NNLO estimate of 
1-T = 0.016±0.003port [24, 26] indicates a much thinner 
tail of the photon energy spectrum and a considerable 
smaller perturbative uncertainty than reported in [20]. 



The DGE analysis thus supports the view that the inte- 
grated photon energy spectrum below i?cut = 1.6 GeV is 
well approximated by a fixed-order perturbative calcula- 
tion, complemented by local OPE power corrections. To 
understand how precisely the tail of the photon energy 
spectrum can be calculated requires nevertheless further 
theoretical investigations. 



IV. NEW PHYSICS IN B ^ Xsj 

Compared with the experimental WA of Eq. (1), the 
new SM prediction of Eq. (4) is lower by more than 1 a. 
Potential beyond SM contributions should now be prefer- 
ably constructive, while models that lead to a suppression 
of the 6 — > S7 amplitude are more severely constrained 
than in the past, where the theoretical determination 
used to be above the experimental one [3] . 

NP affects the initial conditions of the Wilson coeffi- 
cients of the operators in the low-energy effective the- 
ory and might also induce new operators besides those 
already present in the SM. Complete NLO matching 
calculations are available only in the case of the two- 
Higgs-doublet models (THDMs) [27, 28], the minimal su- 
persymmetric SM (MSSM) with minimal-flavor-violation 
(MFV) for small and large tan/3 [28-33], and left-right 
(LR) symmetric models [28]. In the general MSSM [34], 
extra dimensional models like minimal universal extra 
dimensions (mUED) [36] or Randall-Sundrum (RS) sce- 
narios [37], and littlest Higgs (LH) models without [38] 
and with T-parity (LHT) [39] , the accuracy is in general 
strictly LO and hence far from the one achieved in the 
SM. The main features and results of recent analyses of 
beyond SM physics in B ^ Xgj are listed in Tab. I. In 
the following we will briefly review the most important 
findings. 



Even though the effect of charged Higgs boson contri- 
butions in the THDM type II model is necessarily con- 
structive [27, 28], the lower bound on Mh± following 
from B —^ remains in general stronger than all other 
direct and indirect constraints. In particular, B ^ Xgj 
still prevails over B tv [40] for all values of tan /3 apart 
from those lying in the range tan/3 ^ [45,65]. This is il- 
lustrated in the upper panel of Fig. 2. The derived 95% 
confidence level (CL) limit amounts to Mu± > 295 GeV 
independently of tan/3 [13]. In the THDM type I model, 
the strongest constraint on Mj}± stems from the ratio of 
the widths of the Z-boson decay into bottom quarks and 
hadrons, Ri,, and not from B Xg^. 

_ In the MFV MSSM the complete NLO corrections to 
B Xs^ arc also known. The needed two-loop diagrams 
containing gluons and gluinos were evaluated in [28, 29] 
and [30, 31], respectively. Since EW interactions affect 
the quark and squark mass matrices in a different way, 
their alignment is not RG invariant and MFV can only 
be imposed at a certain scale /xmfv that is related to 
the mechanism of supersymmetry (SUSY) breaking [31]. 
For /iMFV much larger than the SUSY masses Msusy, 
the ensuing large logarithms can lead to sizable effects 
in B Xsj, and need to be resummed by solving the 
RG equation of the flavor-changing gluino-quark-squark 
couplings. 

In the limit of A/susy ^ A^u-, SUSY effects can be ab- 
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FIG. 2: Top: Direct and indirect bounds on M^± in the THDIVI 
type II model as a function of tan/3. The colored areas are 
excluded by the constraints at 95% CL. Bottom: 95% CL limits 
on the compactification scale 1/7? in the mUED model as a 
function of the SM central value and total error. The present 
SM result is indicated by the black square. See text for details. 



sorbed into the coupling constants of local operators in an 
effective theory [32, 33]. The Higgs sector of the MSSM is 
modified by these non-decouphng corrections and can dif- 
fer notably from the native THDM type II model. Some 
of the corrections to B ^ Xgj in the effective theory 
are enhanced by tan/3. As a result, they can be sizable, 
of order as tan /? ~ 1 for values of tan/3 ^ 1, and need 
to be resummed if applicable. In the large tan/3 regime 
the relative sign of the chargino contribution is given by 
— sgn(ylt/i). For sgn(At/i) > 0, the chargino and charged 
Higgs contributions interfere hence constructively with 
the SM result and this tends to rule out large positive 
values of the product of the trihnear soft SUSY breaking 
coupling At and the Higgsino parameter fi. 

In the MSSM with generic sources of flavor violation 
a complete NLO analysis is still missing up to date. Ex- 
perimental constraints on generic b —* s flavor violation 
have been studied extensively [34], and radiative inclu- 
sive B-mcson decays play a central role in these analyses. 



In particular, for small and moderate values of tan/3 all 
four mass insertions (<523)ab with A,B = L,R except 
for {S23)rr are determined entirely by S ^ Xg'y. The 
bounds on the mass insertions {S23)ab corresponding to 
tan /3 = 10 are given in Tab. I. For large values of tan/3 
neutral Higgs penguin contributions become important 
and the constraints from both M^M^ ^^'^ Bg-Bs 

mixing surpass the one from B X^j. The effect of the 
precision measurement of the mass difference AMg [35] 
is especially strong in the case of {52^)rl,rr- At large 
tan P the limits on both mass insertions are now imposed 
by the Bg-Bg mixing constraint alone. 

Since Kaluza-Klein (KK) modes in the niUED model 
interfere destructively with the SM 6 — > 57 amplitude 
[36], B -^■ leads to a very powerful bound on the in- 
verse compactification radius of > 600 GeV at 95% 
CL [41]. This exclusion is independent from the Higgs 
mass and therefore stronger than any limit that can be 
derived from EW precision measurements. The 95% CL 
bound on 1 /R as a function of the SM central value and 
error is shown in the lower panel of Fig. 2. In RS mod- 
els, KK modes enhance the BR relative to the SM [37], 
and the bound on the KK masses is in consequence with 
■^KK ^ 2.4 TeV significantly weaker than the constraint 
that derives from EW precision data. 

The contributions to B ^ Xgj from new heavy vec- 
tor bosons, scalars, and quarks appearing in LH models, 
have been studied in [38] for the original model, and in 
[39] for an extension in which an additional Z2 symmetry 
called T-parity is introduced to preserve custodial SU (2) 
symmetry. While in the former case the new contribu- 
tions always lead to an enhancement of B{B Xgj) [38], 
in the latter case also a suppression with respect to the 
SM expectation is possible [39]. As the found LH effects 
in B ^ are generically smaller than the theoretical 
uncertainties in the SM, they essentially do not lead to 
any restriction on the parameter space. 

An alternative avenue to NP analyses of S — > Xgj 
consists in constraining the Wilson coefficients of the op- 
erators in the low-energy effective theory. This model- 
independent approach has been applied combining var- 
ious B- and X-meson decay modes both neglecting 
[42, 43] and including [33, 44] operators that do not 
contribute in the SM. In particular, in the former case, 
merging the information on _B — > Xsj with the one on 
B Xsl~^l~ [45], one can infer that the sign of the 
b ^ sj amplitude is in all probability SM-likc [46]. In 
the case of the Z-penguin amplitude the same conclusion 
can be drawn on the basis of the precision measurements 
of Rh and the other Z ^ bb pseudo observables [43] . 



V. CONCLUSIONS 

The inclusion of NNLO QCD corrections has lead to 
a significant suppression of the renormalization scale de- 
pendences of the B — > Xs7 branching ratio that have 
been the main source of theoretical uncertainty at NLO. 
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The central value of the SM prediction is shifted down- 
ward relative to all previously published NLO results. 
It is now more than Icr below the experimental aver- 
age. This revives the possibility for explorations of new 
physics contributions to rare flavor-changing B-decay 
processes. The dominant theoretical uncertainty in the 
SM is currently due to unknown non-perturbative effects. 
A reduction of this error, together with a calculation of 
the three-loop matrix elements of the current-current op- 
erators and a better understanding of the tail of the pho- 
ton energy spectrum is essential to further increase the 



power of -B ^ Xg^ in the search for new physics. 
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